turn and classical aspects, and therefore, 
it will be useful for implementation of 
quantum computing. 

Recent advances in quantum informa- 
tion theory have inspired an explosion 
of interest in new quantum algorithms 
for solving hard computational prob- 
lems. Three basic “non-classical” proper- 
ties of quantum mechanics — superposi- 
tion, entanglement, and direct 
tensor-product decomposability — were 
main reasons for optimism about capa- 
bilities of quantum computers and 


quantum communications as well as for 
a new approach to cryptography. How- 
ever, one major problem is keeping the 
components of a quantum computer in 
a coherent state, as the slightest interac- 
tion with the external world would cause 
the system to decohere. Another prob- 
lem is measurement: by the laws of quan- 
tum mechanics, a measurement yields a 
random and incomplete answer, and it 
destroys the stored state. 

This proposed reinterpretation of 
quantum formalism opens up new ad- 


vantages of quantum computers: if the 
Madelung equations are implemented 
on a classical scale (using, for instance, 
electrical circuits or optical devices) , all 
the quantum effects important for com- 
putations would be preserved; at the 
same time, the problems associated with 
decoherence and measurement would 
be removed. 

This work was done by Michail Zak of Caltech 
for NASA’s Jet Propulsion Laboratory. For more 
information, contact iaoffice@jpl.nasa.gov. 
NPO-46731 


© Optimal Padding for the Two-Dimensional Fast Fourier Transform 

Appending data to an optimum length decreases computing runtime. 

Goddard Space Right Center, Greenbelt, Maryland 


One-dimensional Fast Fourier Trans- 
form (FFT) operations work fastest on 
grids whose size is divisible by a power of 
two. Because of this, padding grids (that 
are not already sized to a power of two) 
so that their size is the next highest 
power of two can speed up operations. 
While this works well for one-dimen- 
sional grids, it does not work well for 
two-dimensional grids. 

For a two-dimensional grid, there are 
certain pad sizes that work better than 
others. Therefore, the need exists to 
generalize a strategy for determining 
optimal pad sizes. There are three steps 
in the FFT algorithm. The first is to per- 
form a one-dimensional transform on 
each row in the grid. The second step is 
to transpose the resulting matrix. The 
third step is to perform a one-dimen- 
sional transform on each row in the re- 
sulting grid. Steps one and three both 


benefit from padding the row to the 
next highest power of two, but the sec- 
ond step needs a novel approach. 

An algorithm was developed that struck 
a balance between optimizing the grid 
pad size with prime factors that are small 
(which are optimal for one-dimensional 
operations), and with prime factors that 
are large (which are optimal for two-di- 
mensional operations). This algorithm 
optimizes based on average run times, 
and is not fine-tuned for any specific ap- 
plication. It increases the amount of times 
that processor-requested data is found in 
the set-associative processor cache. Cache 
retrievals are 4-10 times faster than con- 
ventional memory retrievals. 

The tested implementation of the algo- 
rithm resulted in faster execution times 
on all platforms tested, but with varying 
sized grids. This is because various com- 
puter architectures process commands 


differently. The test grid was 512x512. 
Using a 540x540 grid on a Pentium V 
processor, the code ran 30 percent faster. 
On a PowerPC, a 256x256 grid worked 
best. A Core2Duo computer preferred ei- 
ther a 1040x1040 (15 percent faster) or a 
1008x1008 (30 percent faster) grid. 

There are many industries that can 
benefit from this algorithm, including 
optics, image-processing, signal-process- 
ing, and engineering applications. 

This work was done by Bruce H. Dean, 
David L. Aronstein, and Jeffery S. Smith of 
Goddard Space Flight Center. Further infor- 
mation is con tained in a TSP (see page 1 ). 

This invention is owned by NASA, and a 
patent application has been filed. Inquiries 
concerning nonexclusive or exclusive license 
for its commercial development should be ad- 
dressed to the Patent Counsel, Goddard Space 
Flight Center, (301 ) 286-7351. Refer to GSC- 
15678-1. 


© Spatial Query for Planetary Data 

This technology is extensible to Earth science and satellite monitoring and surveillance. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Science investigators need to quickly 
and effectively assess past observations of 
specific locations on a planetary surface. 
This innovation involves a location- 
based search technology that was 
adapted and applied to planetary sci- 
ence data to support a spatial query ca- 
pability for mission operations software. 

Conventional databases of planetary 
datasets are indexed and searchable by 
various metadata, such as acquisition time, 


phase of mission, and target. Searching 
these datasets will produce enormous 
datasets that are difficult, or impractical, 
to browse through to identify observations 
of very specific targets. For queries at spe- 
cific locations, it is fundamentally more ef- 
ficient to specify the location as the target 
of the query; and to have the database 
search based on the location of the data 
rather than metadata that is only indirectly 
or tangentially related to location. 


High-performance location-based 
searching requires the use of spatial data 
structures for database organization. Spa- 
tial data structures are designed to organ- 
ize datasets based on their coordinates in 
a way that is optimized for location-based 
retrieval. The particular spatial data struc- 
ture that was adapted for planetary data 
search is the R+ tree. The R+ tree 
arranges data as a set of nodes that repre- 
sents bounding rectangles. Every leaf 


32 


NASA Tech Briefs, April 201 1 


node in die ti'ee is a particular datum 
with coordinates. The root node repre- 
sents a bounding box containing the en- 
tire set of data. Each level of the tree sub- 
divides the search space into smaller and 
smaller bounding rectangles, each con- 
taining a smaller subset of the data. A 
query on the R+ tree for a set of coordi- 
nates will follow one unique path from 
the root to the leaf, and return the data 
contained in the coordinates. The com- 
plexity of the search is bounded by the 
depth of the tree, so the R+ tree insertion 


algorithm maintains a balanced tree such 
that the depth is minimized. 

A map view provides an intuitive way to 
specify a set of coordinates for a location- 
based query. The software will let the 
user select any location on the map of 
the rover’s traverse path to date and re- 
turn all of the results in under a second. 
This was particularly useful during tacti- 
cal activity planning when the question 
of “what data do we have of this location” 
is asked every day. The wrong answer to 
this question is very expensive: false neg- 


atives result in missed science opportuni- 
ties that may never come again, while 
false positives results in a waste of space- 
craft time and bandwidth resources. 

This work was done by Khawaja S. Shams, 
Thomas M. Crockett, Mark W. Powell, Joseph 
C. Joswig, and Jason M. Fox of Caltech for 
NASA’s Jet Propulsion Laboratory. 

The software used in this innovation is 
available for commercial licensing. Please 
contact Daniel Broderick of the California In- 
stitute of Technology at d.anielb@caltech. edu. 
Refer to NPO-46637. 


©Higher Order Mode Coupling in Feed Waveguide of a Planar 
Slot Array Antenna 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A simple technique was developed 
to account for the higher order mode 
coupling between adjacent coupling 
slots in the feed waveguide of a planar 
slot array. The method uses an equa- 
tion relating the slot impedance to 
the slot voltage and a reaction inte- 
gral involving the equivalent mag- 
netic current of the slot aperture and 
the magnetic field coupled from an 
adjacent slot. 

Most waveguide-fed slot antennas use 
centered-inclined coupling slots in the 
feed waveguides. In the proposed 
method, one uses the Elliott’s design 
technique to determine tilt angles and 


lengths of the coupling slots. The radiat- 
ing slots are modeled as shunt admit- 
tances, and the coupling slots are mod- 
eled as series impedances. 

By using reaction integrals to ac- 
count for higher order mode coupling 
between adjacent coupling slots, the 
“active impedance” of each coupling 
slot is obtained, which differs from the 
original impedance by a small amount 
due to higher order mode coupling. 
The tilt angle and length of each cou- 
pling slot are perturbed so that the ac- 
tive impedance of each slot is the same 
as the corresponding starting value of 
the impedance before the iterative 


process started. This process con- 
verges after three or four iterations 
and uses algebraic equations and fast 
reaction integrals. 

Because of desirable features of wave- 
guide arrays such as low volume, ease of 
design and manufacture, and ease of de- 
ployment, these antennas find applica- 
tions in aeronautical and space activities. 
In addition, these all-metal antennas can 
withstand high radiation environments 
encountered in space. 

This work was done by Sembiam Rengara- 
jan of Caltech for NASA’s Jet Propulsion Lab- 
oratory. For more information, contact iaof- 
fice@jpl.nasa.gov. NPO-47419 


@ Evolutionary Computational Methods for Identifying Emergent 
Behavior in Autonomous Systems 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A technique based on Evolutionary 
Computational Methods (ECMs) was de- 
veloped that allows for the automated op- 
timization of complex computationally 
modeled systems, such as autonomous 
systems. The primary technology, which 
enables the ECM to find optimal solu- 
tions in complex search spaces, derives 
from evolutionary algorithms such as the 
genetic algorithm and differential evolu- 
tion. These methods are based on biolog- 
ical processes, particularly genetics, and 
define an iterative process that evolves pa- 
rameter sets into an optimum. 

Evolutionary computation is a method 
that operates on a population of existing 
computational-based engineering mod- 


els (or simulators) and competes them 
using biologically inspired genetic opera- 
tors on large parallel cluster computers. 
The result is the ability to automatically 
find design optimizations and trades, 
and thereby greatly amplify the role of 
the system engineer. 

The ECM technique is extremely ef- 
fective at quickly finding not only indi- 
vidual best solutions, but also an entire 
collection of best solutions for all possi- 
ble ranges of conditions. ECM provides 
an efficient standard of performance 
that specifies the best solution given a 
set of requirements or goals for all possi- 
ble environmental variables. The ECM 
derived standard of performance can be 


directly compared to actual perform- 
ance to provide a test standard of refer- 
ence. Comparison of the ECM best solu- 
tions with actual performance will 
identify where the deviations from opti- 
mum occurred leading to localization of 
problem areas. This includes illuminat- 
ing deficiencies in simulation capabili- 
ties. Examination of the full set of ECM 
derived best performance will enable 
the identification and characterization 
of unexpected or emergent behaviors. 

This work was done by Richard J. Terrile 
and Alexandre Guillaume of Caltech for 
NASA’s Jet Propulsion Laboratory. For more 
information, contact iaoffice@jpl. nasa.gov. 
NPO-46686 
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